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MICROFLUIDIC BULK FLOW DETERMINATIONS 
BASED ON CONVERTING HEAT TRACER MEASUREMENTS 

TECHNICAL FIELD 

[0001] The invention relates generally to methods and systems for 
determining fluid flow and more particularly to measuring flow rates within 
a microfluidic passageway. 

BACKGROUND ART 

[0002] Microfluidic technology is used in systems that perform chemical 
and biological analysis, as well as chemical synthesis, on a much smaller 
scale than previous laboratory equipment and techniques. Microfluidic 
systems offer the advantages of only requiring a small sample of analyte or 
reagent for analysis or synthesis, and dispensing a smaller amount of waste 
materials. A microfluidic system may be a single component or part of a 
larger system. For example, the system may include interface elements 
such as an electrospray ionization tip, which would allow the system to be 
interfaced to a mass spectrometer. The term "microfluidic" as used herein 
refers to features that are fabricated on the micron or submicron scale. For 
example, a typical channel or chamber of a microfluidic system has at least 
one cross-sectional dimension in the range of approximately 0.1 microns 
to 1000 microns. 

[0003] As microfluidic systems increase in complexity, the importance of 
reliable electronic and software processing support to enhance the analysis 
capabilities also increases. Known microfluidic systems provide processing 
support for performing operations such as measuring the flow rate of fluid 
through a system passageway. An accurate determination of flow rate is 
important in a number of different applications, such as in HPLC (High 
Performance Liquid Chromatography) coupled to mass-spectrometer 
analysis. 
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[0004] U.S. Pat. No. 6,386,050 to Yin et al. describes a system and 
method for measuring flow rate within a fluid-bearing passageway of a 
microfluidic system. Heat tracers (i.e., thermal fluctuations) are introduced 
into the flow, so that passage to an interrogation region may be detected and 
timed. The heat tracers may be introduced using an optical heat generator 
or an electrical element, such as a heating resistor. Optical or electrical 
properties of the fluid may be monitored to detect passage of heat tracers 
into the interrogation region. The Yin et al. patent teaches that the flow rate 
is based directly upon calculating the speed of a heat tracer. 

[0005] A number of other systems for measuring flow rates are described 
in the Yin et al. patent. For example, in the described system of U.S. Pat. No. 
4,938,079 to Goldberg, intrusive monitoring is used. In the Goldberg system, 
the heat tracer is introduced into the flow of liquid by a microwave heating 
device or a source of focused infrared energy. The flow rate is measured by 
determining the transit time of the heat tracer from its source to a sensor. 
As one possibility, the dielectric constant of the liquid is monitored to detect 
changes in temperature. The Goldberg system is designed to provide 
accuracy at flow rates of less than 100 cc/hour. U.S. Pat. No. 5,726,357 
to Manaka and U.S. Pat. No. 5,623,097 to Horiguchi et al. also describe 
microfabricated devices which employ thermal approaches to calculating 
flow rates. 

[0006] While the known approaches to calculating flow rates based on 
measuring the time interval required for a heat tracer to travel a particular 
distance provide useful information, more accurate calculations are sought. 

SUMMARY OF THE INVENTION 

[0007] In a microfabricated device, the on-chip fluid flow rate is calculated 
by incorporating an adjustment that is specific to the scaling between a bulk 
fluid flow rate and a propagation rate of a heat tracer introduced into a 
particular microfluidic passageway. Thus, a chip-based microfluidic device 
may include a heat generator coupled to introduce the heat tracer into a flow 
of fluid through the passageway, a detector positioned to detect temperature- 
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dependent variations of the fluid within an interrogation region of the passage- 
way, and a processor that is configured to determine the transit time of the 
heat tracer through the passageway and to determine the bulk fluid flow rate 
on the basis of the transit time and a pre-identified scaling between heat 
tracer propagation rates and bulk fluid flow rates of the microfluidic passage- 
way. 

[0008] In accordance with the method, the heat tracer is introduced into 
the fluid flowing through the microfluidic passageway of the microfabricated 
device. A temperature-dependent condition of the fluid is monitored during 
passage through the interrogation region. On the basis of the monitoring, the 
arrival of the heat tracer within the interrogation region is detected. The bulk 
flow rate of the fluid through the passageway is then calculated using the 
heat tracer transit time and the pre-identified scaling between the propagation 
rate of the heat tracer and the bulk fluid flow rate. In order to enhance the 
accuracy and/or the dynamic range of flow rate measurements, a series of 
microfluidic features (e.g., sensing electrodes) may be formed along the 
microfluidic passageway at increasing distances from the source of the heat 
tracer, so that repetitive measurements can be acquired for a single heat 
tracer. 

[0009] The scaling in the two rates is a function of the properties of the 
fluid and the properties of the structure directly related to forming the 
microfluidic passageway. The relevant properties include the heat capacities 
of the fluid and the structure, as well as the geometry of the structure which 
forms the passageway. The propagation rate of the heat tracer will be less 
than the flow rate of the fluid. The calculation of the appropriate scaling factor 
may include determining the product of the ratio of the two heat capacities 
times the ratio of a cross-sectional dimension of the microfluidic passageway 
to the cross-sectional dimension of the structure that forms the passageway. 
For example, the second ratio may be the half width of the passageway to the 
half width of the passageway-forming structure. 
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[001 0] In application, the heat capacity of the passageway-forming 
structure and the relevant dimensions of the structure are known and fixed. 
Therefore, only the heat capacity of the fluid varies, where the heat capacity 
is a function of the specific heat and the density of the fluid. In one embodi- 
ment, the conversion from the propagation rate information of the heat tracer 
to the bulk fluid flow rate is determined algorithmically. As an alternative, the 
conversion is performed by applying a calibration curve which implements an 
upward adjustment to the bulk fluid flow rate. 

[001 1] The heat generator for introducing the heat tracers and the detector 
for monitoring the temperature-dependent variations in the fluid may be inte- 
grated onto the same substrate on which the passageway and other micro- 
fluidic features are microfabricated. The microfluidic features may include 
micromachined components, such as valves. Moreover, the processor and 
the storage of the passageway data and fluid data necessary for the conver- 
sion may be integrated onto the same substrate as the microfluidic features. 
Tracer generation may be optically based or electrically based. Similarly, 
tracer detection may be performed electrically or optically. 

BRIEF DESCRIPTION OF THE DRAWINGS 

[0012] Fig. 1 is a perspective view of an example of a microfabricated 
system having a flow rate conversion capability in accordance with the 
invention. 

[001 3] Fig. 2 is a schematic representation of on-chip components of the 
system of Fig. 1 . 

[0014] Fig. 3 is a schematic representation of one possible implementation 
of a heat generator for introducing a heat tracer into a microfluidic 
passageway of the system of Fig. 1 . 

[001 5] Figs. 4 and 5 are graphs of the conductivity of fluids as a function of 
temperature. 
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[001 6] Fig. 6 illustrates transit times of heat tracers at different flow rates. 

[001 7] Fig. 7 is an illustration comparing the expected bulk fluid flow rate 
to the actual bulk fluid flow rate. 

[001 8] Fig. 8 is a sectional view of a microfluidic passageway of Fig. 1 . 

[001 9] Fig. 9 is a process flow of steps for carrying out flow rate conver- 
sions in accordance with one embodiment of the invention. 

DETAILED DESCRIPTION 

[0020] With reference to Fig. 1 , a microfabricated device 10 is shown as 
being simplified to a single microfluidic feature, namely a microfluidic pas- 
sageway 12 having an inlet 14 and an outlet 16. In a more complex embodi- 
ment, the substrate 18 in which the passageway 12 is fabricated includes an 
interconnection of microfluidic features, such as transfer channels, separation 
channels, mixing channels, and microfluidic compartments which may be 
used for fluid treatment and analysis. The fluid treatments may include, but 
are not limited to, mixing, labeling, filtering, extracting, precipitating, digesting, 
and the like. The substrate 18 may also include microfluidic flow control 
elements such as mixers, valves, pumps, pressure regulators and mass flow 
regulators. 

[0021] The microfabricated device 10 may be formed by direct means, 
such as photolithographic processing, wet or dry chemical etching, laser 
ablation, or traditional machining. The device 10 may also be fabricated by 
indirect means, such as injection molding, hot embossing, casting, or other 
processes that utilize a mold or patterned tool to form the microfluidic 
features. 

[0022] The substrate 18 may be formed of a polymer, glass, silicon or 
ceramic material. Polymers offer advantages, with polyimide being particu- 
larly well suited. A polymer may be selected from the following classes, 
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polyimide, PMMA, polycarbonate, polystyrene, polyester, polyamide, poly- 
ether, polyolefin, and mixtures thereof. 

[0023] As previously defined, the term "microfluidic" refers to a component 
or feature associated with channels or chambers that are generally fabricated 
on the micron or submicron scale. For example, the microfluidic passageway 
12 of Fig. 1 may have a cross-sectional dimension in the range of 0.1 microns 
to approximately 500 microns. While the inlet 14 and the outlet 16 create the 
impression that the passageway is circular, conventional fabrication tech- 
niques result in non-circular passageways. 

[0024] The microfluidic passageway 12 of Fig. 1 is shown as dashed lines, 
since the passageway is formed in the interior of the multi-layer substrate 18. 
At an upstream side of the passageway is a pair of electrodes 20 and 22 on 
opposite sides of the passageway. Referring briefly to Fig. 3, the electrodes 
may be thin metallic films that are patterned during the fabrication of the 
microfluidic system. The electrodes are connected to a source 24 of alternat- 
ing current by traces 26 and 28. While the source 24 is shown as applying a 
sine wave, square wave activation or other current fluctuations may be sub- 
stituted. 

[0025] In Fig. 2, the mechanism for introducing the heat tracers into the 
microfluidic passageway 12 is identified broadly as a "heat generator" 30. 
The heat generator may include a controller that receives control signals from 
a processor 32, such as a microprocessor bonded to the substrate 18. The 
controller may include transistors or other electronic components which are 
integrated with other elements of the microfluidic system during fabrication. 
That is, the electronic components may be formed using photolithography 
or other integrated circuit fabrication techniques. However, a less complex 
embodiment utilizes bonding techniques, such as flip-chip bonding, to con- 
nect the processor 32 and any other electronic components to the system. 
As another possibility, the operations performed by the electronic components 
may be provided off-chip, so that only electrical connections to the substrate 
18 are required. 
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[0026] As an alternative to the use of the electrodes 20 and 22 and the AC 
source 24 of Fig. 3, the heat generator 30 may be a serpentine pattern of 
traces through which direct current is conducted to generate heat in the area 
adjacent to the microfluidic passageway 12. Microfabricated heaters are well 
known in the art. While there are benefits to providing contactless conduction 
of heat to the fluid within the passageway, the heat generator 30 may use an 
approach in which the heated element, such as a heating resistor, is in direct 
contact with the fluid. 

[0027] The thermal energy that is introduced to the fluid within the micro- 
fluidic passageway 12 by the heat generator 30 should be sufficiently great 
that the change in temperature of the fluid is perceivable using techniques to 
be described below. However, the increase in temperature should not 
degrade the chemical composition of the fluid. Optionally, the system may 
include a feedback loop to the controller of the heat generator 30 in order to 
provide active temperature control. The same concerns are applicable if the 
electrical heating is replaced with optical heating, such as the use of a lamp 
or laser having a wavelength that is absorbed by the fluid within the passage- 
way. 

[0028] As previously noted, the processor 32 may be a prefabricated 
component that is bonded to the substrate 18. The processor contains 
software and/or firmware that, through its operations, guides or controls the 
actions of the integrated circuitry. For example, the processor may contain 
programmable logic that allows a programmed algorithm to be executed so 
as to perform certain functions, including the conversion from a heat tracer 
propagation rate to a bulk fluid flow rate. Bonding the processor can involve 
utilizing contact solder to connect corresponding electrical contact points on 
the processor and the substrate 18. The contact solder may be tailored to the 
maximum temperature that will be withstood by the microfluidic processing. 
While not shown, the processor or the substrate may include contact pads 
that are used to connect the system to remote (i.e., off-system) devices or 
systems. 
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[0029] The substrate 18 may also include micromachined components, 
such as the valve 34. Techniques for micromachining valves and other 
mechanical devices onto a substrate are known, particularly if the substrate 
is formed of a semiconductor material. 

[0030] A detector 36 is aligned with an interrogation region 38 of the 
microfluidic passageway 12. Prior to activation of the heat generator 30, the 
fluid within the passageway will have a particular level of conductivity. The 
conductivity will depend upon the ionic species that are constituent to the 
fluid, upon the electrolyte that is used, and upon the temperature of the fluid. 
The detector monitors the condition of the fluid flowing through the interroga- 
tion region, so that changes that are induced as a consequence of activations 
of the heat generator are detected. This approach can be used in any appli- 
cation in which there is a movement of the fluid (as compared to applications 
in which separations of constituents do not require bulk fluid movement, e.g., 
electrophoresis). The approach is particularly well suited for chromatographic 
separations and to the accurate dispensing and general manipulation of 
liquids. 

[0031] In some applications, the microfabricated system is used for on-chip 
analysis, for example, separations based on sample hydrophobicity, charge, 
chirality, or affinity to a specific compound. In other applications, additional 
analysis may be performed following the release of the fluid from the sub- 
strate 18, such as schematically represented by the electrospray 40 in the 
direction of a charge plate 42. As one example, the on-chip components may 
implement an HPLC separation of the sample, and the electrospray may 
soft-ionize the sample for introduction into a mass spectrometer. 

[0032] In an optical embodiment of the detector 36, a mirror may be 
placed on a side of the microfluidic passageway 12 opposite to a non-thermal 
optical source such as a laser. Prior to heating of the fluid by the heat 
generator 30, the light from the non-thermal laser will undergo a degree of 
refraction upon entering and exiting the passageway. Thus, the axis of the 
light that reaches the detector will have a first angle relative to the axis of the 
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light as it approaches the passageway. The detector may be a position- 
sensitive detector, such as a charged coupled device (CCD). On the other 
hand, when a heat tracer reaches the interrogation region 38 of the passage- 
way, the refractive index will change and light that reaches the CCD will have 
a second angle that is different from the normal state (i.e., unheated state) of 
the fluid. 

[0033] Rather than monitoring the refractive index of the fluid in order to 
determine when a heat tracer reaches the interrogation region 38, the 
conductivity of the fluid can be monitored. Fig. 1 shows a two-electrode 
embodiment of the conductivity monitoring. In this embodiment, electrodes 44 
and 46 are fabricated on opposite sides of the microfluidic passageway 12. 
Four-electrode embodiments may also be used, as described in U.S. Pat. No. 
6,386,050 to Yin et al. A voltage generator, such as an AC source, is con- 
nected to two of the electrodes. The frequency of the voltage generator may 
be one megahertz, but other frequencies may be used. The generator 
provides signals that are 180 degrees out of phase. The electrodes may be 
conductive traces on the substrate 18. At least one capacitive cell is formed. 
Since the electrical conductivity of the fluid within the passageway will affect 
the capacitive cell, the conductivity can be monitored by appropriately con- 
nected sense circuitry. The electrical conductivity of the fluid will vary with 
temperature, so that arrival of a heat tracer at the interrogation region is 
detectable. 

[0034] As opposed to the planar electrodes 44 and 46 represented in 
Fig. 1, the electrodes may be rings that are formed at the interrogation region. 
The rings are spaced apart and the electrical conductivity of the fluid between 
the rings is monitored. Similarly, co-planar electrodes that are spaced apart 
along the direction of fluid flow may be used to monitor electrical conductivity 
within the region between the electrodes. 

[0035] While only a single interrogation region is shown in Figs. 1 and 2, 
more reliable flow rate information may be possible using multiple interroga- 
tion regions. For example, rather than determining a heat tracer transit time 
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as the time between activating the heat generator 30 and detecting the 
temperature-dependent fluctuation within the interrogation region 38, transit 
time may be measured as the time required for a previously introduced heat 
tracer to pass from an upstream interrogation region to a downstream inter- 
rogation region. Accuracy may be further enhanced by forming a series of 
electrodes (or other sensing features) at increasing distances from the source 
of the heat tracers. As a consequence, repetitive measurements may be 
acquired for each heat tracer, allowing a greater dynamic range in flow rate 
measuring as compared to using a single electrode at a fixed position from 
the source of the heat tracer. 

[0036] In practice, the distance over which the heat tracer travels in 
determining the transit time may be 10 millimeters (mm). A heat tracer may 
be introduced by activating a heater with a heat pulse of between one and 
two volts DC at a current between 0.35 and 0.7 amperes. The approach of 
detecting temperature-dependent variations within the interrogation region 38 
may be conductivity based, rather than optically based, such as by detecting 
in-phase current using a lock-in amplifier. Fig. 4 is a plot 48 of the electrical 
conductivity of a solution of 0.1 percent Formic Acid in 50 percent Acetonitrile 
(CH 3 CN) as a function of temperature. The conductivity is represented in 
standardized units (uS/cm), with the conductivity increasing as the tempera- 
ture rises from 14 degrees Celsius to 41 degrees Celsius. Similarly, Fig. 5 is 
a plot 50 of the electrical conductivity of a solution of 0.1 percent Formic Acid 
in water. Again, conductivity increases with temperature. 

[0037] On-chip mobile phase linear velocity measurements were 
acquired or a solution having a mobile phase of 0.1 percent Formic Acid 
and 25 percent Acetonitrile. Fig. 6 shows information for determining four 
different transit times for four pump flow rates, which are measured in 
microliters per minute (uL/min). In each of four instances 52, 54, 56 and 58, 
the application of a heat pulse is marked by a baseline offset 60 at the 
initiation time. The Y axis is not set to scale in Fig. 6. Conductivity is 
monitored within the interrogation region and remains substantially at the 
baseline until a baseline change 62, 64, 66 and 68 occurs. The conventional 
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(prior art) technique is to calculate the on-chip fluid flow rate as a direct 
correlation to the transit time of the heat tracer. 

[0038] In Fig. 7, a graph 70 of the expected velocity value as a function of 
flow rate is shown. However, Fig. 7 also includes four plotted values of 
measured flow rate velocity as a function of the pump flow rate. The four 
measured values are connected by a linear fit line 72, so as to better illustrate 
the difference between expected velocities and actual velocities. It has been 
determined that there is a linear scaling between the propagation rate of a 
heat tracer and the bulk fluid flow rate. The cause of the lag of the heat tracer 
propagation rate is represented in Fig. 8. A thermally generated heat tracer 
74 within a fluid 76 exchanges thermal energy with the walls 78 and 80 that 
form the microfluidic passageway 12. The exchange of thermal energy is 
represented by the arrows 82. As a consequence of the exchange, the walls 
78 and 80 will include heated regions 84 and 86. While the bulk fluid flow 
rate is not affected, the flow rate of the heat tracer 74 will be reduced by the 
thermal exchange. The degree of effect is a function of the heat capacities 
of the fluid and the walls, as well as the geometries of the passageway and 
walls. Specifically, the propagation rate of the heat tracer 74 is a function of 



1 + — Pw °p w 
a Pf ° P f 

where p w is the density of the wall, p f is the density of the fluid, c pw is the 
specific heat of the wall, c pf is the specific heat of the fluid, a is the inside 
radius (or half width) of the microfluidic passageway 12, and w is the outside 
radius (or half width) of the structure that forms the microfluidic passageway. 
This equation is generally accurate for two-dimensional considerations. 
Scaling factors for a three-dimensional application can be calculated through 
numerical calculations for the heat transfer. 

[0039] As an example of the use of the equation, f will be equal to approxi- 
mately 0.35 if w = 50 fjm, a = 10 fjm, the fluid is water, and the material for 



Attorney Docket No. 10030605-1 



12 

forming the passageway is Kapton, which is a registered trademark of E.I. Du 
Pont for a flexible film of electrically insulative material. Table 1 includes 
some of the information relevant to this determination. 



TABLE 1 


Property 


water 


Kapton 


p (kg/m 3 ) 


997 


1420 


k (W/m-K) 


0.608 


0.12 


C p (J/kg-K) 


4177 


1090 



[0040] As can be seen, f reflects the percentage of the heat tracer 
propagation rate to the bulk fluid flow rate. The percentage is a function of a 
geometry ratio (i.e., w/a) and a heat capacity ratio. 

[0041] Fig. 9 is a process flow of steps for executing the invention. At 
step 88, information relevant to determining flow rates for the microfluidic 
passageway is stored. Thus, the wall thickness, the passageway dimensions, 
and the specific materials that form the walls are relevant. Thus, a table may 
be stored so as to include the different scaling factors for a given structure 
and for a variety of fluids that might be used. Alternatively, an equation can 
be stored, where the equation results from a fit to numerical calculations that 
are performed as a function of the heat capacity of a fluid within the structure. 
As one possibility, the storage occurs within memory of the processor 32 of 
Fig. 1. Alternatively, the substrate 18 may have separate storage capacity, 
such as integrated read only memory (ROM). As another possibility, the 
information is stored off-chip for access when the conversion from heat tracer 
propagation information to bulk fluid flow information is implemented. The 
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conversion may occur off-chip, but there are advantages to enabling all 
processing to take place at the chip level. 

[0042] In some applications, there will only be a limited number of fluids 
that will be used in the analysis. For example, there may be only two fluids 
contemplated for use within the system of Fig. 1 . Optionally, the structure that 
forms the microfluidic passageway 12 may be selected to be relatively neutral 
with respect to the selection of the fluids in a particular use. That is, the 
structure may be designed to minimize the heat partitioning effect such that 
the sensitivity to the type of fluid is low. This may simplify the processing. 

[0043] At step 90, a heat tracer is introduced into the flow of fluid through 
the microfluidic passageway 12. Heat pulses may be applied to the narrow 
zone along the passageway through which fluid is flowing. A heat tracer may 
be introduced via a resistive heater element placed along the passageway or 
by means of propagating an electrical current directly through a narrow zone 
of fluid between two electrodes 20 and 22. Optical heating is also an option. 

[0044] As the heated volume of fluid flows through the microfluidic pas- 
sageway 12, the heat partitions into and out of the walls 78 and 80 which 
define the microfluidic passageway 12, as described with reference to Fig. 8. 
As a consequence, the heat tracer propagates at a slower rate than the bulk 
fluid flow rate. As indicated at step 92, the interrogation region is monitored 
to detect arrival of the heat tracer. The interrogation region is monitored to 
detect temperature-dependent variations of the fluid. Optical monitoring may 
occur by measuring the refractive index of the fluid. On the other hand, 
electrical conductivity may be monitored using either direct contact electrodes 
or "contactless" detection, as described in U.S. Pat. No. 6,386,050 to Yin et al. 

[0045] At steps 94 and 96, the temperature-dependent variation within the 
interrogation region is detected and the transit time of the heat tracer over a 
known distance is determined. The known distance may be the distance 
between a heating zone and the interrogation region or between a pair of 
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interrogation regions. Then, the conversion from the heat tracer propagation 
information to the bulk fluid flow information is performed at step 98. 

[0046] The conversion to the bulk fluid flow rate may be performed directly 
from the transit time determined within step 96. As an alternative, the heat 
tracer propagation rate may be calculated and then used in the conversion 
process. The previously identified algorithm provides the basis for converting 
from the slower heat tracer propagation rate to the bulk fluid flow rate. As 
compared to an algorithmic conversion, the information stored at step 88 may 
include a calibration curve for converting measured heat tracer time intervals 
to known bulk fluid flow rates. 

[0047] An advantage of the processing of Fig. 9 is that the flow rate dif- 
ferential has been recognized and a methodology has been incorporated to 
accurately determine the bulk fluid flow rate. Accurate calculations of flow 
rates are significant to some analytical processes. 



